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The electronic structure of the vacancy-ordered Ko.5Fei.75Se2 iron-selenide compound (278 phase) 
is studied using the first-principles density functional method. The ground state of the 278 phase is 
stripe-like antiferromagnetic, and its bare electron susceptibility shows a large peak around (tt,tt) 
in the folded Brillouin zone. Near Fermi level, the density of states are dominated by the Fe-3d 
orbitals, and both electron-like and hole-like Fermi surfaces appear in the Brillouin zone. Unfolded 
band structure shows limited similarities to a hole doped 122 phase. With O.le electron doping, 
the susceptibility peak is quickly suppressed and broadened; while the two-dimensionality of the 
electron- like Fermi surfaces are greatly enhanced, resulting in a better nesting behavior. Our study 
should be relevant to the recently reported superconducting phase Ko.5+xFei.75+ySe2 with both x 
and y very tiny. 

PACS numbers: 75.25.-j,71.20.-b,71.18.+y 



Introduction- The discovery of vacancy-ordered iron- 
selenides 1 4 has greatly renowned the study of iron-based 
superconductors 5 . These compounds share a general 
chemical formula of A y Fe2- x Se2, and their crystal struc- 
tures are close to BaFe2As2- However, they differ from 
the latter because the iron sites are not fully occupied in 
the selenides, and the iron- vacancies may form ordered 
structure at certain x. For example, in the vacancy- 
ordered K0.8Fe1.6Se2 (or K2Fe4Se5, 245 phase), the iron- 
vacancies form a superstructure of yfb x y/E (in unit of 
nearest neighbor iron-iron distance), separating square 
blocks of four iron- atoms. The distinct geometry and the 
magnetic frustration thereafter lead to a unique antiferro- 
mangetic (AFM) ground state pattern, i.e. the block-spin 
AFM24 

The presence of the ordered vacancy patterns alters the 
symmetry of the crystal lattice, and therefore its effect 
is far beyond the simple rigid-band shifting model. The 
new crystal symmetry modulates the electronic structure 
of the vacancy- free AFe2Se2 (or 122) phase, and the iron 
vacancies introduce extra scattering centers that further 
modify the electronic states. Therefore, the electronic 
structure of these vacancy-ordered phases can be com- 
pletely different one from another, although their chem- 
ical composition may look fairly close. It therefore leads 
to a question that among all the possible vacancy-ordered 
phases, which one is the actual parent phase for the su- 
perconductivity. As the 245 phase seems to be present 
universally in these superconducting samples, it was pro- 
posed that the pairing was mediated by AFM fluctuation 
in the doped 245 compound^. This argument was chal- 
lenged by a phase-separation scenario that the 245 phase 
and the superconducting phase are spatially separated. 
The phase separation argument has been supported by 
angle resolved photo-emission spectroscope (ARPES)&2 
and scanning tunneling microscope (STM) or transmis- 
sion electron microscopy (TEM) experiments^ - — , which 
suggest the vacancy concentration is much smaller and 



are disordered in the superconducting phased. 

Very recently, there has been an experimental report 
on T C =33K superconducting iron selenides, which con- 
sists of superconducting bricks embedded in a back- 
ground of 245 phase. The bricks have chemical com- 
ponent Ko.5+ccFei.75 +2y Se2. In these bricks which take 
up about 20% of the volume, it is found that x ^0.15 
and y ^0.01. This raises an interesting possibility 
that the vacancy-ordered Ko.5Fei.75Se2 (278 phase here- 
after), may be a parent compound of the superconducting 
phase 14 . The authors presented the preliminary verifica- 
tion on this possible parent phase by conducting self- 
consistent measurements of chemical compositions and 
the atomic imaging through scanning tunneling. There- 
fore, it will be interesting and urgent to provide first 
principles calculations and theoretical studies on this 278 
phase. Note that unlike the 245 phase, the single crys- 
tal 278 phase has not been discovered yet, and the phase 
form only "spider-web" like filaments. 

In this article, we present our latest first principles 
results of the 278 phase. We examine the structural dis- 
tortions induced by the iron vacancies; calculate its band 
structure and density of states, and compare them with 
the 122 phase. We present and compare its bare elec- 
tronic susceptibility and nesting functions at zero and 
O.le doping levels. From these calculations, we conclude 
that the ground state for the 278 phase is stripe-like 
AFM and metallic. Upon O.le electron doping, which 
quickly suppresses the long range AFM order, the two- 
dimensionality of the Fermi surface sheets are greatly en- 
hanced, resulting in a better nesting behavior. 

Method- In the current study, we employed den- 
sity functional theory (DFT) using plane-wave basis 
as implemented in Vienna Abinit Simulation Pack- 
age (VASP) 15 i 16 . The ion- valence electron interactions 
were modeled using projected augmented wave (PAW) 
method; and the Perdew, Burke, and Ernzerhof flavor of 
exchange-correlation functional were chosen^. A 480 eV 
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energy cut-off for the plane- wave basis and a 6x6x3 T- 
centered K-grid ensures the total energy converges less 
than 1 meV/atom. The lattice constants and inter- 
nal atomic positions were fully optimized until forces on 
individual atoms are smaller than 1 meV/A and total 
stress less than 0.05 kB. The density of states (DOS) 
calculations were performed on a much denser K-grid of 
24x24x12 and tetrahedra method. 

The obtained band structure were then fitted to a 
tight-binding hamiltonian of 70 Fe-3d orbitals using max- 
imally projected wannier function method^. The result- 
ing hamiltonian was then used to interpolate the elec- 
tronic states on a dense T-centered 100 x 100 x 100 K-grid. 
The bare electronic susceptibility xo is then calculated 
using 
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where e^k and /(e M k) are the band energy (measured 
from Ep) and occupation number of |/ik), respectively; 
\n) denotes the n th wannier orbital; and iVk is the number 
of the k points used for the irreducible Brillouin zone 
(IBZ) integration. It is worthy noting that one should not 
make the usual assumption that (fik\n) = 1, as pointed 
out by S. Graser et a/.—. The information is also used to 
calculate the nesting function using 
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Results and discussion- Figure [T(a)1 and [T(b) 
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unit cell we chose to perform our calculations. While the 
iron vacancy order is crucial to the electronic structure 
of these compounds, the effect of potassium order seems 
to be negligible except for the resulting chemical poten- 
tial shifting, as shown in previous studies 3,4 . Therefore, 
without losing generality, we chose a unit cell with high 
symmetry where occupied pot assiu m sites are closest to 
the iron- vacancy sites. Figure |l(c) compares the IBZ of 
a 122 phase unit cell (enclosed by blue dotted lines) and 
the one of a 278 phase unit cell (enclosed by red dotted 
lines). It is worthy noting that the area of the first Bril- 
louin zone (BZ) for the 122 phase (the blue square) is 
four times the one for the 278 phase (the red hexagon). 

The optimized lattice constants for the 278 phase are 
a = 8.3813 A and c = 13.9846 A (Fig. [TJbl, respectively. 



Previous density functional studies of iron-pnictides and 
iron-selenides severely suffer from a c-collapse problem 
that optimized lattice constant c are significantly smaller 
(beyond usual DFT limit of 5%) than the experimental 
value. This particular problem is in fact absent in the 278 
phase if one assumes that the 278 phase c value should 
not be too far from the 122 phase (14.0367 A). The 7 an- 
gle formed by a and b lattice vectors is 127.12°, which is 
slightly larger than 126.87°. The latter is the ideal value 
of 7 if the iron square lattices are perfectly preserved. 
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FIG. 1: (a) The unit cell of 278 phase. The red atoms are 
iron; the blue atoms are potassium; and the yellow atoms are 
selenium, (b) The iron-plane showing the ordered vacancies 
and the choice of lattice vectors as in (a). It is obvious that the 
two lattice vectors are both vTO (in unit of nearest neighbor 
Fe-Fe distance), (c) The first BZ (solid lines) and IBZ (dotted 
lines) of the 122 phase (the blue lines) and the 278 phase (the 
red lines). Both M' and X' in the 122 IBZ are X in the 278 
IBZ. 



Therefore, the square lattice of iron is slightly distorted 
by less than 0.5%. Furthermore, the optimized lattice 
constants can be roughly translated into a' = 3. 748 A in 
the 122 phase, which is 4.3% smaller than the experimen- 
tal value, suggesting a contraction of lattice constants 
due to the vacancy ordering. The presence of the ordered 
vacancies also introduces further distortions of the atomic 
positions within one unit cell. The seven iron-sites in one 
unit cell splits into three nonequivalent classes: the four 
iron atoms that are nearest neighbors of the iron- vacancy 
(Fe 1 ); the iron site that is in between two neighboring 
iron- vacancies (Fe 11 ); and the rest two iron atoms (Fe m ), 
and their fractional heights are 0.2488, 0.2463 and 0.2523, 
respectively. In addition, the Fe 1 atoms slightly move to- 
wards the nearest vacancy site while both Fen and Fem 
remains in square lattice position (so that they only dis- 
place in z-direction). The Fe-Se-Fe angles formed by the 
next-nearest-neighbor iron-sites and the above selenium 
atom are in the range of 108.16° to 109.66°, slightly larger 
than the one in the KFe 2 Se 2 122 phase (106.61°). 

With the optimized crystal s truct ure, we th en study 
its electronic structure. Figure |2(a) and 2(b) show the 
band structure and density of states (DOS) of 278 phase. 
As the unit cell of the 278 phase is 8 times the 122 phase 
primitive cell (which is not orthogonal), its irreducible 
Brillouin zone (IBZ) as well as the band structure is 
heavily folded (for the IBZ folding, please refer to Fig. 
l(c)[ ). However, it can still be identified that three bands 
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FIG. 2: (a) Band structure and (b) density of states of 278 
phase. Unlike other iron-pnictides or chalcogenides, the elec- 
tronic states of 278 phase from Ef-3.0 eV to Ef+2-0 eV are 
almost exclusively from Fe-3d orbitals. 



FIG. 3: Band structure of (a) 278 phase in the unfolded BZ 
(122 phase BZ), and (b) 122 phase in 122 phase BZ. In both 
figures, the Fermi level Ef is aligned at zero energy (red line), 
while in the 122 phase plot, an additional -0.25 eV line (cor- 
responding to the doping level in the 278 phase) is shown. 



crosses the Fermi level, suggesting three Fermi surface 
sheets in the BZ. Unlike either 122 or 245 phase, the Fe- 
3d bands in the 278 phase are well separated from all 
other contributions, and dominates the electronic states 
from E^-3.0 eV to £"^+2.0 eV as s uggested in the pro- 
jected DOS calculations (Fig. |2(b)[ ). Several flat bands 
can be identified close to Ef-, from T to M and from Z to 
A, which contributes the DOS peak at Ef- It is worthy 
noting that the flat band contribution to conductivity is 
almost negligible due to the large effective mass. Thus, 
the pseudogap-like feature from Ep-1-0 eV to Ef +1.0 
eV may be consistent with the "bad- metal" behavior ob- 
served in the experiments. 

As the electronic states near Ef are almost exclu- 
sively contributed by the Fe-3d orbitals, we can easily 
fit its band structure within £p-3.0 eV to Ef+2.0 eV to 
a tight-binding hamiltonian using maximally projected 
wannier orbitals. The original 278 phase band structure 
is then unfolded to the 122 phase IBZ using the method 
introduced by W. Ku et al^. The resulting unfolded 
band structure is shown in Fig. |3(a)| in comparison with 
the 122 phase band structure (Fig. |3(b)[ ). Due to the 
strong structural distortion and iron vacancies, the band 



structure of the 278 phase is significantly smeared out, 
although some feature of the latter can still be identified 
from the unfolded 278 band structure. In particular, both 
the upper and lower band limit of Fe-3d orbitals seem to 
be shifted « 0.25 eV, meaning the Fermi level Ep of the 
278 phase is shifted « —0.25 eV. It is worthy noting that 
this shifting of Ep corresponds to « 0.5 \e\ hole doping 
in the 122 unit cell, or Ko.75Fe2Se2 in chemical formula. 
One can also identify some other characteristic rigid-band 
shifting in the unfolded 278 phase band structure, espe- 
cially along Y-Z since the direction is mostly unaffected 
by the iron- vacancies whose effect is primarily in-plane. 

Features beyond rigid-band shifting can also be identi- 
fied in the unfolded band structure. First of all, the addi- 
tional scattering centers caused significant band smear- 
ing and splitting in the regions around Ef — 1.0 eV and 
Ef + 1-0 eV. Secondly, The -0.25 eV Fermi level shifting 
allows two hole-like Fermi surface sheets emerge around 
T in the 122 phase, as one can clearly identify two hole 
like band dispersion cross the -0.25 eV line in Fig. [3(b)] 
This feature, however, is absent in the 278 phase band 
structure. One another shallow hole-like band crossing, 
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however, can be identified around X in the 278 phase 
unfolded band structure, which is clearly missing in the 
rigid-band shifted 122-phase. Finally, the bands around 
Ep are extremely smeared out except from T to Z, sug- 
gesting that the electron transport suffers from strong 
scattering in the a-b plane and may exhibit a low con- 
ductivity. Due to the complexity of the band structure, 
it is rather difficult to conclude the type of the major 
carrier, but both electron and hole seems to play roles. 
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FIG. 5: Fermi surface sheets of (a) parent 278 phase, and (b) 
O.le electron-doped 278 phase. The F point is at the center 
of each plot. 
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FIG. 4: (a-b) Bare electron susceptibility xo and (c-d) Nesting 
function / of parent 278 phase (a, c) and O.le electron doped 
278 phase (b,d). The plots are made in a full BZ with F at 
four corners. 

The bare electronic susceptibility xo is also calculated 
using the maximally projected wannier function hamilto- 
nian. Fig. |4(a) and |4(b)| show xo at k z = plane for both 
the parent 278 phase and O.le-doped phase, respectively. 
Previous calculations for BaFe2As2 and LaOFeA o 2Q i 21 
show significant enhancement at (jr, 0) in the two-iron 
(unfolded) BZ. For the parent 278 phase, a similar en- 
hancement of ^100% can also be identified around (tt, tt) 
in the folded BZ, suggesting an AFM long range order. 
In fact, mean- field calculations based on the extended 
J i — J 2 model have proposed four possible ground state 
candidates 23 . By employing spin-polarized density func- 
tional calculations, we confirm that the ground state par- 
ent 278 phase is stripe-like order AFM, which is ^150 
meV/Fe lower than the non-magnetic phase and ~30 
meV/Fe lower than the second lowest AFM order. With 
O.le electron doping, the enhancement is greatly reduced 
and broadened, thus the long range AFM order may be 
destroyed at such doping level. 

Finally, we present the Fermi surfaces for both the par- 
ent compound and the O.le doped compound (Fig. [5]). 
As pointed out previously, despite of the complex band 
structure of the 278 phase, only three bands across the 



Fermi level, forming three Fermi surface sheets. Upon 
electron doping, the sizes of the two hole pockets along 
T-M are dramatically reduced while the electron pock- 
ets are greatly enlarged. Apart from that, a fourth sheet 
emerges around F . Furthermore, the two-dimensionality 
of both electron pockets are significantly enhanced, lead- 
ing to a better nesting behavior. This is also reflected in 
the nesting function plots (Fig. 4(c)| and |4(d)[ ), which 
show over twice enhancement of the nesting function 
around the center of BZ. 

Conclusion In conclusion, we have performed first- 
principles study on vacancy-ordered Ko.sFei.75Se2 (278 
phase). The ground state of the parent 278 phase is 
stripe-like AFM and metallic, whose magnetic order may 
be destroyed by electron/hole doping. The Fe-3d or- 
bitals dominate the electron DOS near Fermi level, and 
both electron- and hole-like Fermi surfaces emerge in 
the BZ. The unfolded band structure shows limited fea- 
tures of a hole-doped 122 phase, but effects beyond rigid- 
band model must be taken into consideration to under- 
stand the electronic structure of the 278 phase. Fi- 
nally, a O.le electron doping significantly enhances the 
two-dimensionality of its Fermi surfaces, resulting in a 
better nesting behavior. In our present calculations, 
electron correlations beyond local density approximation 
level have been neglected. A strong on-site Coulomb re- 
pulsion could alter the magnetic order and lead to Mott 
insulator, in light of the AFM Mott insulator in 245 com- 
pound. The analysis is beyond the scope of the present 
paper and experiment synthesis of the 278 phase would 
shed new light on this issue. 
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